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Abstract For a e]0, 7r/2[. let (S)q be the control system i = (F + uG)x, where x belongs to the 
U: two-dimensional unit sphere S^, u £ [—1,1] and F, G are 3x3 skew-symmetric matrices generating 

■ rotations with perpendicular axes of respective length cos(a) and sin(Q;). In this paper, we study 

the time optimal synthesis (TOS) from the north pole (0, 0, 1)"^ associated to as the parameter 

a tends to zero. We first prove that the TOS is characterized by a "two-snakes" configuration on 
^ ' the whole S*^, except for a neighborhood Ua of the south pole (0, 0, —1)^ of diameter at most 0{a). 

We next show that, inside Ua, the TOS depends on the relationship between r{a) := 7r/2a— [7r/2a] 
and a. More precisely, we characterize three main relationships, by considering sequences (Q!A:)fc>o 
satisfying (a) r{ak) = r; (b) r{ak) = Cak and (c) r{ak) = 0, where f £ (0, 1) and C > 0. In each 
^ I case, we describe the TOS and provide, after a suitable rescaling, the limiting behavior, as a tends 

. to zero, of the corresponding TOS inside Ua- 
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1 Introduction 

Let a g]0, 7r/2[. On the unit sphere S'^ C M'^, consider the control system (S)q defined by 

{T,)a x={F + uG)x, X ^ {xi,X2,xzf , = 1, |m| < 1, (1) 

where F and G are two 3x3 skew-symmetric matrices representing two orthogonal rotation s wi th axes of 
length respectively cos(a) and sin(a), a €]0, 7r/2[ (for the precise meaning of length, see Section 2^). With no 
loss of generality, we assume that 



G := " sin(a) . (2) 






In this paper, we aim at describing the time optimal synthesis (TOS for short) from the north pole N := 
(0,0,1)"^ for i.e. for every x G we want to find the time optimal trajectory steering TV to a; in 

minimum time (see Figure 0) . 

In particular we are interested in the qualitative shape of the time optimal synthesis in a neighborhood of 
the south pole S = (0, 0, —1)"^, in the limit a — *■ 0. The interest for that problem stems from motion planning 
issues in aeronautics and quantum control, see Q for instance. 

The present paper is actually a continuation of ||^ in the sense that it answers questions raised in the latter 
paper. There, the purpose was to provide a lower and an upper bound for N(a), the maximum number of 
switchings for time optimal trajectories for the left invariant control system 

(5)„ g = giF + uG), <? G 50(3), |m| < 1, (3) 

where F and G are defined in (||). Recall that, for such control systems, it is known (cf. for instance ^ |^) 
that every time optimal trajectory is a finite concatenation of bang arcs (i.e. u = ±1) or singular arcs (u = 0). 
A bang arc is an integral trajectory corresponding to the rotations 

X+:= F + G, X^ :=F- G, (4) 

and is denoted by e'^^^a;, t € [0,T], where e = ±, a; is the starting point of the bang arc and T is its time 
duration. Moreover, a switching time - or simply a switching - along a time optimal trajectory is a time so 
that the control u is not constant in any open neighborhood of to- 
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To estimate N{a), a suitable Hopf map 11 : SO{3) S*^ was introduced to project {S)a onto (S)q. In 
particular, every time optimal trajectory of is the projection by 11 of a time optimal trajectory of {S)a- 
It results that, if a time optimal trajectory on has a certain number of switchings, then this number is lower 
than or equal to the maximum number of switchings for the optimal problem on SO{3). The construction of 
time optimal trajectories of (S)q, was performed according to the general theory of time optimal synthesis on 
2-D manifolds developed in §, |, |, |, |l^, |ll|, |l|l , and recently gathered in the book 

The question of studying N{a) was first addressed in ||l| where, using the index theory developed by 
Agrachev, the authors proved that N{a) < [Tr/a], where [•] stands for the integer part. That result was not 
only an indirect indication that N{a) would tend to infinity as a tends to zero, but it also provided a hint on the 
asymptotic of N{a) as a tends to zero. Notice that for a = the systems (]l|) and (^) are not controllable. With 
the techniques developed in enough properties for the TOS associated to a < 7r/4, were identified in 

order to improve the upper bound of [nj and to actually show that, for a small 



TV (a) < kn j + 5, where k 



M 
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In it is proved that, for a < n/A, the extremals associated to (S)q, (i.e. the trajectories candidate for time 
optimality obtained after using the Pontryagin Maximum Principle -PMP for short-), starting from the north 
pole N are bang-bang trajectories, i.e. finite concatenations of bang arcs of the type 

, where the initial time duration Si verifies Si £ (0, tt], all the time durations of the interior bang arcs are equal 
to v{si), where the function v is defined in Eq. (|l2|) below, and the final time duration Sf verifies Sf < v{si). 
Of particular importance for the construction of the TOS, are the switching curves, i.e. the curves made by 
points where the control switches from +1 to —1 or viceversa and defined inductively by 

dis) = e^-^("'e^-'"iV, Cl{s) = e^^^'^'^C^^^is), (where e = ±1 and fc 2, fcA/). (5) 

Since the PMP gives just a necessary condition for optimality, it is crucial to determine the time after which 
an extremal is no more optimal. In |^, we showed that the number of bangs must be lower than or equal 
to fcj\/ + 1 and the extremals cover the sphere according to the "two-snakes" configuration as depicted in 
Figure The two "snakes" correspond to extremal trajectories starting respectively with control +1 and — 1. 
For more details, see 

However, in 1^], we were not able to construct the complete TOS associated to (E)^. In particular, we could 
not show the optimality of all the extremals up to kM — 1 bangs arcs and we could not complete analytically the 
construction of the synthesis in a neighborhood of the south pole S. There, the minimum time front develops 
singularities due to the compactness of S'^. We only provided numerical simulations describing the evolution 
of the extremal front in a neighborhood of the south pole. As a ^ 0, these numerical simulations suggested 
the emergence of an interesting phenomenon (see Fig. ^ : define the remainder 

r{a) := 7r/2a- [7r/2a]. (6) 

Then, there are three possible patterns of TOS in the neighborhood of the south pole S, each of them 



depending on a relation between r{a) and a. See Section 2.2 below, where these relations are formulated as 
well as conjectures. 

In §, the TOS for (S„) was studied in the context of quantum control. The control system (Sq) describes 
the population transfer problem for the xs-component of the spin of a (spin 1/2) particle, driven by a magnetic 
field, which is constant along the xa-axis and controlled along the .Ti-axis, with bounded amplitude. In that 
paper, the TOS, for a > tt/4 was completed and, in the case a < 7r/4, further information was obtained, 
for what concerns time optimal trajectories steering the north to the south pole (in fact the most interesting 
trajectories for the quantum mechanical problem). Such optimal trajectories belong to a set S containing at 
most 8 trajectories, half of them starting with control -f-l and the other half starting with control —1, and 
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F-G F+G 




Figure 2: The "two-snakes" configuration defined by the extremal flow. Notice that this set of trajectories 
covers the whole sphere, but in principle not all extremals are optimal and a point can be reached by more 
than one trajectory at the same or at different times. 



switching exactly at the same times. It was also proved that the cardinality of S depends on the remainder 
r(a) defined in Eq. (^). For instance, for a and r(a) small enough, then S contains exactly 8 trajectories (four 
of them are optimal) while if r{a) is close to 1, then 5 contains only 4 trajectories (two of them are optimal). 

The purpose of the present paper consists in studying the TOS associated to (S)q as a tends to zero, 
focusing in particular on its behavior inside a neighborhood of the south pole. Roughly speaking, we want to 
determine, as a tends to zero, what could be a possible limit for the TOS associated to (S)^ (as suggested 
for instance by the patterns depicted in Fig. ^) and then to prove the convergence (in some suitable sense) of 
the TOS associated to (S)q to that limit. To proceed, we embark on the study of a geometric object T{a,T) 
called the extremal front at time T along (S)q and defined as the set of points reached at time T by extremal 



trajectories starting from N (see Section 3.1 for a precise definition). The extremal front J-{a,T) contains the 
minimum time front OF(a,T), i.e. the set of points reached at time T by time optimal trajectories. When 
J^{a,T) = OF{a,T), we say that !F{a,T) is optimal. 

We first prove, in the case in which fcjv/ is odd (being the other case analogous), that the extremal front 
J^{a, kMT^) is made up of the union of two curves E^{a, •) : (0, tt] 5^, e = ±, with £^(a, •) = Ilajgf ""^(a, •), 
where I^x^ is the orthogonal symmetry with respect to the xs-axis. Moreover, for a small enough, £^(q;, •) 
admits a convergent power series of the type X]/>o fii^i'^'i'^))'^^ i where the f[{s,r) are real-analytic functions 
of (s,r) G E^, 27r-periodic in s with 

/0\ /-2rcA / f(4r + c,)s2 \ 

/o+(^, ''^ " I °J ' ^1^^"' '^^ " I ^''^ I ' ^2''^'' ^) = I + ^2.).s. 1 . (7) 

As a trivial consequence, we deduce that for r G [0, 1], s G R and a small enough, we have 

S%a, s) - /o^(,s, r{a)) + fl{s,r{a))a + f'^{s,r{a))a' + 0{a^). (8) 

and 

^£^(a,5) = ^Jl{sA<^))^ + §;r2{sAc^)W +0{a^). (9) 
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Figure 3: Conjectured shapes of the synthesis in a neighborhood of the south pole. Switching curves are 
curves made by points in which the control switches from +1 to —1 or viceversa. Overlap curves are curves 
made by points reached optimally by more than one trajectory. The curve 7fc is a bang arc that is also an 
overlap curve since trajectories having a different history travel on it at the same time. 
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where [©(a'^)] < Clap with C >0 constant mdepcndent of {r,s,a). 

Then we show that !F{a,T) is actuaUy optimal for T < {knj — l)7r and a small enough (see Remark 7 
below). Moreover, we show that ^(a, (fcjv/ — l)7r) is a circle of radius 2{l + r{a))a up to order (see Remark 3 
below). As a consequence of the optimality of J^{a, (fcA/ — l)7r), we get that all the extremals of the "two-snakes" 
configuration depicted in Figure ^ are optimal up to time (fcjvf — l)'"'). In other words, if Ua is the connected 
component of S*^ \ J-^{a, {kM — l)7r) containing the south pole, we obtain the optimal synthesis on S'^ \Ua. 
Notice that Ua is a neighborhood of the south pole of size proportional to a. 

The expressions (^-(§|) are central tools to understand the possible asymptotic behaviors of the TOS 
associated to (E)a, as a tends to zero. 

For this purpose we observe that the expressions of and f2 in (^, depend explicitly on the remainder 
r{a). This fact suggests the need to impose particular relationships between a and r{a) in order to define 
any asymptotic behavior. In other words we must let a goes to zero, only along certain subsequences {ak)k>o 
where a specific relationship holds between at and r(afe). The analysis of Eq. (|^) will help us to determine 
such relationships and to prove that the conjectures made in |^] about the qualitative shape of the synthesis 



near the south pole were true (see Section 2.2 and Figure In particular we will see that there are exactly 
three qualitatively different asymptotic behaviours of the synthesis as a goes to zero, described by the following 
cases. 

First, we analyze the case in which a is arbitrarily small, with r(a) G (0, 1) uniformly far from and 1. To 
simplify further the discussion, it is reasonable to consider the following. 

(CI) For f G (0, 1), let a tend to zero along the subsequence ak ■— 2{k+f) ' ^'^ f{ak) = r. 

In this case £'^(a, •) is approximated, up to order a^, by the expression S+ff{-, f). As a consequence J-{a, kMT^) 
is approximately a circle of radius 2fa centered at the south pole. We are then able to give a qualitative 
description of the optimal synthesis, as stated below in Theorem This synthesis turns out to be exactly the 
one described in Figure || (case B), as predicted in Q. 

It remains then to study the cases in which r(ct) can be arbitrarily close to or 1. For this purpose we first 
consider the case in which r{a)la remains bounded above and below by positive constants as a tends to zero. 
From Eq. (|^) it is clear that this is equivalent to say that /f(-,r)a is comparable to f2{-,r)a'^. For simplicity 
we consider the following. 

(C2) For C > 0, let a tend to zero along a subsequence (Q;fc)fe>o such that r{ak) = Cafe. 

In this case £^{a, ■) is well approximated by + (/f (■, C) + /|(-, 0))a^. If C > 7r/4, the synthesis is equivalent 
the one of the previous case. On the other hand if C < 7r/4 the synthesis is more complicated (sec Section ||) 
and it turns out to be exactly the one described in Figure ^ (case C), as predicted in j^. 

If a and r(a) tend to zero with r(a)/a tending to infinity (resp. to zero) it is possible to see that the 
synthesis is qualitatively equivalent to the one of case (CI) (resp. (C2)). 

The third interesting case is the following. 

(C3) Let a tend to zero along the subsequence ak ■= so that r(afc) = 0. 

In this case the extremal front at time kMT^ contains the south pole and the corresponding optimal front reduces 
to that point. The optimal synthesis is then described starting from the extremal front ^(a, (^a/ — l)7r) = 
OF {a, (kM — l)""), and it corresponds to the one described in Figure ^ (case A), as predicted in p. 

With similar arguments, one can see that in the case in which a is small and r{a) is close to 1, the optimal 
synthesis is qualitatively equivalent either to that of Case (CI) or to that of Case (C3), and this concludes the 
description of the possible asymptotic behaviors as a tends to 0. 

Remark 1 It is interesting to notice that numerical simulations show that for a decreasing to zero continuously, 
the qualitative shape of the optimal synthesis described in Figure || alternates cyclically in the order BCABCA.... 

Let us describe the results obtained in the case (CI) in more details. Since J- {a, kf,[-K) is approximated, 
up to 0{a^), by a circle of center S and radius 2fa, we are able to show that it is optimal, so that all the 
extremals of the "two-snakes" configuration depicted in Figure Hare optimal up to time fcj\/7r. In other words. 
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Figure 4: Optimal synthesis for the linear pendulum 



if Va is the connected component of \ J- {a, k^ii^) containing the south pole, we obtain the optimal synthesis 
on 5*2 

As a tends to zero, Va collapses on S. Hence one must rescale the problem by a factor 1/a, in order 
to describe the TOS inside Vq. Also notice that since we are in a neighborhood of the south pole we can 
project the problem on the plane (xi, X2)- We are now in a position to define a possible limit behavior for the 
TOS inside Va- Let Ma be the linear mapping from M.^ onto defined as the_composition of the projection 
(xi, a;2, Xa) ^ {xi,X2) followed by the dilation by \/a. Denote by (S)q (resp. OF{a,kMi^)) the image by Ma 
of (resp. OF(a, kMT^))- Then, {T,)a is a perturbation by 0{a^) of the standard linearized pendulum 

(^en) : ( = ""^2' , . , , ^ (10) 

^ ^ \ Z2 = zi+u, (21,2:2) SK , |w| < 1 ^ ^ 

while OF{a, kMT^) is a perturbation by 0{a^) of C(0, 2f), the planar circle of center (0, 0) and radius 2f. As a 
consequence, the candidate limit TOS inside Va is the one associated to the problem of reaching in minimum 
time every point of the ball B(0, 2f) starting from C(0,2f), along the dynamics of the standard linearized 
pendulum. To prove such a result, we first study the above mentioned optimal control problem and show 
that the corresponding TOS is characterized by an overlap curve 7°g„, which is the set of points z e with 
212:2 > and belonging to the locus (see Figure 0) 

zf + zl + 2zlzl - if^zf + (4 - 4f2)2| = . 

The optimal synthesis inside C(0, 2f) is then described by the following feedback, defined on B{0,2f) \ 7° : 
"above" 7pe„, the control u is constantly equal to — 1 a nd "below" 7°g„, it is constantly equal to 1 (see Fig. ^). 



Finally, the asymptotic result we prove in Section 4.2 is the following 



Theorem 1 For r E (0,1), let {ak)k>i be the sequence defined by := 2{k+f} f'-'^ k > 1. Consider 7°g„, 
the overlap curve of the TOS for the optimal control problem consisting of starting from C(0, 2f), the planar 
circle of center (0,0) and radius 2f, and reaching in minimum time every point of B(0,2f) along the control 
system (jJ^. Then, for k large enough, the TOS associated to {'S)ak inside OF{ak,kM'!^) is characterized by 
an overlap curve 7°^ so that the optimal feedback takes the value —1 "above" 7°^, and the value 1 "below" 
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. Moreover, 7°^ converges to 7°g„ in the topology, uniformly with respect to r in any compact interval 
0/(0,1). 

The resuhs in the cases (C2) and (C3) are described in more details in Sections ^ and |[ 

Remark 2 Notice that the sequence {ak)k>i defined above has been chosen in order to simpUfy the previous 
statement. Indeed the same result could be restated in a more general way by taking an arbitrary sequence 
(Q!fc)fe>i converging to zero and such that r{ak) converges to f, or letting the remainder vary on a compact 
subinterval of (0, 1). 

The paper is organized as follows. In the second section, we collect basic facts, notations, results, and 
conjectures of Q. The third section gathers the detailed description of the extremal front and the proof of 
Eq. (|^). Sections ^ || and ^ treat respectively the cases (CI), (C2) and (C3). In the appendix, we finally 
prove a technical result needed throughout the paper. 



2 Notations and previous results 
2.1 Basic Facts 

Definition 1 An admissible control u{.) for the system is a measurable function u{.) : [a, 6] [—1, 1], 

while an admissible trajectory is an absolutely continuous function x{.) : [a,b] — )■ S'^ satisfying a.c. for 
some admissible control u{.). If x{.) is an admissible trajectory and u{.) the corresponding control, we say that 
{x{.),u{.y) is an admissible pair. 

For every x G S'^, the minimization problem consists of determining an admissible pair steering the north pole 
to X in minimum time. More precisely 

Problem (P) Consider the control system For every x G S'^, find an admissible pair (x(.),u(.)) 

defined on [0,T] such that x{0) ~ N, x{T) ~ x and x{.) is time optimal. 

An optimal synthesis from the north pole (in the following optimal synthesis, for short) is the collection of 
all the solutions to the problem (P). More precisely 

Definition 2 (Optimal Synthesis) An optimal synthesis for the problem (P) is the collection of all time 
optimal trajectories T = {xx{.) ■ [0, 62] 1-^ 5*^, x G S"^ : Xx{0) = N, Xx{bx) = x}. 

For more elaborated definitions of optimal synthesis see |l^ and references therein. The standard tool to look 
for optimal trajectories is a first order necessary condition for optimality known as the Pontryagin Maximum 
Principle (PMP for short), (cf. |l3)) as stated below for our minimum time problem on S^. 

Define the following real- valued map on T*S'^ x [—1,1], called Hamiltonian, 

n{X, X, u) :=< A, {F + uG)x > . 



Set: 



H{X,x) := max H{X,x,v). (11) 

vel-is] 



The PMP asserts that, if 7 : [a, b] — > 5*^ is a time optimal trajectory corresponding to a control u : [a, 6] 
[—1, 1], then there exists a nontrivial Held of covectors along 7, that is a never vanishing absolutely continuous 
function X : t G [a,b] 1-^ X{t) G T*^^^S'^ and a constant Aq < such that, for a.c. t G Dom{j), we have: 

i) m = -^(A(t), 7(t), u{t)) - -A(t)(F + uit)G), 

ii) n{X{t),j{t),u{t)) +Xo = 0, 
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iii) n{x{t),j{t),u{t))^H{j{t),x{t)). 



In the more general case in which the target and the initial datum (also called source) are two smooth manifolds 
Ao and A/i the previous statement must be modified by adding the so-called transversality conditions: 

iv) < A(a),iJ >= Vw e T^(<j)7Vo , <A(fe),ui>=0 Vw e T^(6)7Vi. 

Remark 3 A trajectory 7 (resp. a couple (7, A)) satisfying the conditions given by the PMP is said to be an 
extremal (resp. an extremal pair). An extremal corresponding to Ag = is said to be an abnormal extremal, 
otherwise we call it a normal extremal. 

Definition 3 (bang, singular for the problem (0)-(|2|)) A control u{.) : [a,b] [—1,1] is said to be 
a bang control if u{t) = +1 a.e. in [a,b] or u{t) = —1 a.e. in [a,b]. A control u{.) : [a, 6] [—1,1] *s 
said to be a singular control if u{t) — 0, a.e. in [a,b]. A finite concatenation of bang controls is called a 
bang-bang control. A switching time of u{.) is a time t G [a,b] such that, for every e > 0, u is not bang or 
singular on {t — e,t + e)n [a,b]. A trajectory of the control system is said to be a bang trajectory (or 

arc), singular trajectory (or arc), bang-bang trajectory, if it corresponds respectively to a bang control, singular 
control, bang-bang control. If i is a switching time, the corresponding point on the trajectory x{t) is called a 
switching point. 



2.2 Description of Previous Results 

In 1^, it was proved that, for every couple of points, there exists a time optimal trajectory joining them. 
Moreover it was proved that every time optimal trajectory is a finite concatenation of bang and singular 
trajectories. More precisely we have: 

Proposition 1 For the minimum time problem associated to for each pair of points p and q belonging 

to S'^, there exists a time optimal trajectory joining p to q. Moreover every time optimal trajectory for ([|j-(^ 
is a finite concatenation of bang and singular trajectories. 

Notice that the previous proposition docs not apply if a = or a = tt/2, since in these cases the controllability 
property is lost. 

In 1^ it has been proved that a = 7r/4 is a bifurcation for the qualitative shape of the time optimal synthesis, 
for instance the time optimal synthesis contains a singular arc if and only if a > 7r/4. Since in this paper we 
arc interested in the limit a — > 0, in the following wc always assume a < tt /A. In this case, using the PMP, the 
following properties characterizing the optimal trajectories were established in |^: 

i) x{.) is bang bang; 

ii) the duration Si of the first bang arc satisfies Si G (0, tt], 

iii) the time duration between two consecutive switchings is the same for all interior bang arcs (i.e. excluding 

the first and the last bang) and it is equal to v{si), where v[-) is the following function 



sin(s) 
cos(s) + cot^(a) 



v[s) = TT + 2 arctan ( , ) ■ (12) 



One can immediately check that this function satisfies u(0) = v(k) = tt and v{s) > tt for every s e (0, tt) 

iv) the time duration of the last arc is Sf e (0,f(si)], 

Moreover, thanks to the analysis given in one easily gets (always in the case a < 7r/4): 

v) the number of switchings Nj. of a;(.) satisfies the following inequality 



TT 

2^ 



(13) 
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Figure 6: Locally optimal switching curves and non locally optimal switching curves with the corresponding 
synthesis 



Conditions i)-v) define a set of candidate optimal trajectories. Notice that conditions i)-v) are just necessary 
conditions for optimality and one is faced with the problem of selecting, among them, those that are really 
optimal. In particular, given a trajectory satisfying conditions i)-v), one would like to find the time after which 
it is no more optimal. 

Some questions remained unsolved, in particular questions relative to local optimality of the switching curves 
defined in Eq. (^. Roughly speaking we say that a switching curve is locally optimal if it never "reflects" the 
trajectories (see Fig. ^ A).^ When a family of trajectories is reflected by a switching curve then local optimality 
is lost and some cut locus appear in the optimal synthesis. 

Definition 4 A cut locus is a set of points reached at the same time by two (or more) time optimal trajectories. 
A subset of a cut locus that is a connected manifold is called overlap curve. 

An example showing how a "reflection" on a switching curve generates a cut locus is portrayed in Fig. ^ B 
and C. More details are given later. More precisely the following questions were formulated in 0: 



Question 1 Are the switching curves Cf(s), s £ (0, tt], locally optimal? More precisely, one would like to 
understand how the candidate optimal trajectories described above lose optimality. 

Question 2 What is the shape of the optimal synthesis in a neighborhood of the south pole? 

Numerical simulations suggested some conjectures regarding the above questions. More precisely, in Q the 
following conjectures were made: 

^ More precisely consider a smooth switching curve C between two smooth vector field Yi and Y2 on a smooth two dimensional 
manifold. Let C{s) be a smooth parametrization of C. We say that C is locally optimal if, for every s G Dom{C), we have 
C{s) 7^ aiYi{C(s)) + ce2Y2{C'{s)), for every 01, a2 s.t. Jiia2 > 0. The points of a switching curve on which this relation is not 
satisfied are usually called "conjugate points". See Fig. BI 
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CI The curves C^(s), (fc = 1, km) arc locally optimal if and only if fc < [^^j^t] ^ 

Analyzing the evolution of the minimum time wave front in a neighborhood of the south-pole, it is reasonable 
to conjecture that: 

C2 The shape of the optimal synthesis in a neighborhood of the south pole depends on the remainder r{a) 
defined in Eq. (||). Notice that r{a) belongs to the interval [0, 1). More precisely, it was conjectured in 
H] that for a £ (0,7r/4), there exist two positive numbers ai and a2 such that < ai < a < q;2 < 2a 
and: 

CASE A: r{a) S (^, 1). The switching curve C|^^ glues to an overlap curve that passes through the 
origin (Fig. |3|, Case A). 

CASE B: r{a) e The switching curve Cf^^ is not reached by optimal trajectories in the interval 

]0, tt]. At the point C^^^ (0), an overlap curve starts and passes through the origin (Fig. H, Case B) . 

CASE C: r(a) S (0, ^). The situation is more complicated and it is depicted in the bottom of Fig. ^, 
Case C. 

For r = 0, the situation is the same as in CASE A, but for the switching curve starting at C^^^_^{0). 

As explained in the introduction, the presence of several cyclically alternating patterns of optimal synthesis, 
each of them depending on an arithmetic property of a, was already confirmed in [Q, by counting the number 
of optimal trajectories reaching the south pole. 

Remark 4 The first conjecture is implicitly disproved by the results of this paper. More precisely an immediate 
consequence of our results is that the switching curve C|^^_2 is always locally optimal, while C^^^_^ is not, in 
general. However, for every fixed r < ^ there exists a small enough with f < r(a) < ^ such that Ci:^^_i is 
locally optimal too, which contradicts the conjecture. On the other hand Conjecture C2 is correct and, at the 
light of our main results, is completely proved and clarified. 

2.3 Notations 

All along the paper we use the notation e = ±1. The set so(3) of 3 x 3 skew-symmetric matrices is a three- 
dimensional vector space on which the following bilinear map 

< A,B >= -Tr{AB), A,B G so(3), 

is an inner product. For A S so(3), ||^|| := V< A,A> is the norm (or length) of A. With the above notations, 
F and G are perpendicular and normahzed so that = cos(a) and ||G|| = sin(a). 

Let Id be the 3x3 identity matrix. We recall that N ~ (0, 0, 1)"^ and denote the south pole as S* = 
(0, 0,-1)^. Set ct := cos(t) and St := sin(t) for t e [0, 2tt). Recall that X+:= F + G and := F -G and we 
have 

/O -c„ \ /O -c„ 0\ 

X+ = Cq -Sa , = \ Ca Sa\ ■ 

\0 J \0 J 

Let be the orthogonal symmetry with respect to the xa-axis, i.e. is represented in the canonical basis 
by Diag{^l, —1, 1). Then, we have the following trivial but useful property. 

Ux,X, = X-ellx,. (14) 

We next recall standard formulas for a rotation e*^ of SO (3) in terms of its axis Y (whose length is equal 
to one) and its angle t. We have 

e^^ = Id + stY +{l~Ct)Y^. (15) 

Moreover, for t G [0, 27r), we have 

^e(t)z-(t) ._ ^tx+^tx^ ^0(t)z+it) ^tx^^tx+ 
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where the unit vectors Z^(t), Z^(t) arc defined by 




-C{t) 

Z+it) = \ C{t) I , = ( C{t) I , (16) 

^B{t) 

with B{t) := /"^^ ^ , C(0 - , /'^^ , and the angle Q{t) by 

e(t) = 2 arccos(sj/2C2Q - cl^^). (17) 



3 The Extremal Front 
3.1 Definition and description 

As said in the introduction. T{a,T) the extremal front along (E)q at time T is the set of points reached at 
time T by extremal trajectories starting from N, i.e. 

T{a, T) := {x e S'^ : 3 an extremal pair (x(.), A(.)) such that x(0) = iV, x{T) = x}. (18) 

Such extremals are parametrized by the length of the first bang arc, the one of the last bang arc and the number 
of arcs: 

n terms 



S-(s,t) = e^.'*e^^^'"^") ■ ■■e^+"<")e^-", iV, (20) 
n' terms 

where s G (0, tt], t e (0, i'(s)], the number of bang arcs (n and n' respectively) is an integer and 
(-) e = +1 (resp. e = —1), if n is odd (resp. even), 
(-) e' = +1 (resp. e' = —1), if n' is even (resp. odd). 

Roughly speaking, we would like to compute the limit, as a — > 0, of J- (a, T), when T is such that the extremal 
front reaches a neighborhood of the south pole. 

The idea is that, once one knows the extremal front J-{a,T) and if it is optimal, then one can continue to 
build the synthesis for times bigger than T using J- (a, T) as a source for the minimization problem. 

The identification of the front J-{a, T) is not easy since it requires the computation of the product of several 
exponentials of matrices. Moreover, if !F{a, T) crosses some switching curve, then the number of exponentials 
in general depend on the point. 

This problem is overcome by considering J-{a, T) only at times equal to multiples of tt. Indeed, first notice 
that, for T = TT [^], the extremal front reaches the points C^^^(O), i.e. the points where the last switching 
curves C^^^ start. Thanks to Proposition || below, at these times, every extremal trajectory has the same 
number of switchings. The extremal front at times that are not multiple of tt can be obtained a posteriori, 
continuing the extremal front, as explained above. 

From the structure of the extremal trajectories it follows that the time at which the point C^(s) is reached 
is Tk{s) = s + kv{s). 



Lemma 1 Let k be an integer satisfying 1 < k < Mrr 
function of s. 



[cotjaf ~ If 
2col{af - 1 



, then Tk (s) is a strictly increasing 
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Proof. It holds 

J l + 2cs cot{af + cot{af + k (2 + 2c, cot(a)^) 

-^Tk{s) = 5 ^ 2 ^. (21) 

ds l + 2csC0t{af +cot{af 



It is clear that the denominator of the above fraction is never vanishing on [0,7r] if a < 7r/4. On the other 
hand the numerator, as a function of s, reaches its minimum at s = tt, where it is equal to (cot(a)^ — 1)^ — 
fc(2cot(Q!) — 1), and then the conclusion follows easily. ■ 
As a consequence, we obtain the following important corollary. 

Corollary 1 Let k be an integer satisfying 1 < k < Afmon- If o-n extremal trajectory is switching at time 
T — kn, then the length s of the first bang arc satisfies s — n. 

Since for a small /cm < A/'mon, then, for T = kn where A: is a positive integer such that k < [7r/(2a)], we 
have that all the extremal trajectories switch exactly k times (except the trajectories with length of the first 
switching equal to tt that switch fc — 1 times). Therefore, the extremal front J- {a, kn) is described by the next 
proposition. 

Proposition 2 Let k be a positive integer such that 1 < fc < [71/(20;)]. Then, if a is small enough, we have 

T{a,kTT) = {£+{a,k,s), s e (0, tt] } |J {£" (a, fc, s), se(0,7r]}, where: (22) 

r g(fc.^(fc-iM.)-.)x_gii^ieM.))z_(.(.))g.x+^ fork odd, , . 

1 ^{k(^-v{s))-s)X+^^e{v{s))Z-{v{s))^sX+]^ fork even. ^ ' 

The expression for £^ is the same as the expression for £^ after exchanging the subscripts + and —. As a 
consequence, = n^,^^, where IVx^ is the orthogonal symmetry with respect to the x^-axis. 

Remark 5 Notice that £'^{a, fc, 0) — £^^{a, fc, tt), e = ± , so that is described by a continuous closed curve. 

3.2 Description of the extremal front jF(a, /ca/ti") and consequences 

As sketched in the introduction, we must describe the optimal synthesis on deprived of a neighborhood of 
the south pole. For that purpose, we will provide the precise asymptotics of J^{a,kMT^), as a tends to zero, 
and derive, from its topological nature, the minimum time front at time fcA/7r. 

From now on, for simplicity, we drop the dependence of £^ on kM, i.e. we set £^ {a, s) := £^(a, fcM, s), and 
we assume that fcj\/ is odd. 

In the following, it will be useful to think of a and r as two independent variables. For this purpose, define 

tp{a, r, s) := (5^ - r) (tt - v{s)) + v{s) - s , 

0ia,r,s) := (JL _ i±l) 9(^,(5)) , 

X^(a, r, s) := e.^i°'-'^''')^-^Q^i°''^-'^)^-^i^i^))e^^^ TV . 

It is clear from (^) that 

£'{a,s) = x^(a,f(a),s) . 
The following result is the key point in order to describe the extremal front at time fcA/7r. 

Lemma 2 There exists ao > such that the function , e = ±, defined above, is real-analytic for (r, s, a) G 
X /, where I = (— ao,ao)- Moreover, it admits a convergent power series 

X^(a,r,s)-^/f(s,r)a', (24) 

i>0 
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where the ff{s,r) are real- analytic functions of {s,r) S K^, 2'K-periodic in s (therefore they are bounded over 
Rx [0,1];. 

As a consequence, the extremal front J^(a, /cmtt), which is a continuous closed curve, is piecewise analytic 
with discontinuities at s = 0,tt for derivatives of order greater than or equal to one. 

Proof of Lemma |2|. 

We will prove the proposition only for ■ Since is 27r-periodic in s and r enters in an afRne way in ip 
and 0, the real issue of analyticity revolves around the variable a. First of all, it is clear that w(s) is actually 
a real-analytic function for (s, a) G M x /, where / = (— ao, cto) with ao > small enough. Therefore, one has 
only to prove the real-analyticity of rp{a, s) := and where (3{s, a) :— Q(v{s)), for (s, a) G M x /, 

where / = {—ao, ap), for some ao > 0. 

Note that ^ 

s) = — arctan(s^/i(s)), /i(s) 



a " c2 + s^Cs 

The function /i is real-analytic for (s, a) G Kx/, with / open neighborhood of zero, and thus uniformly bounded 
over M X /. In addition, arctan(-) is real analytic in a neighborhood of zero. Hence the conclusion for 'il'{a, s). 
As for rewrite first Eq. @ as 

cos(/3(s)) = 1 - G(s,a), 

with 



G{s,a) := 2si 



(25) 



We first need to determine a convergent power series for (3 from the expression 

/3 = arccos(l - G). (26) 

Note that |G(s, a)| < 5a^ for a small enough. We first expand arccos(l — G) in a power series in G. Starting 
from the power series 

(l-t)-'/' = l+ 

m>l 

with radius of convergence equal to 1 we get 

-^(arccos(l — G)) = ^= — , ^ , 

dG^ ^ " V2G^l-G/2 

and, after simple integration, 

arccos(l - G) = - V2G(l + V ttG™) , (27) 

m>l ^ ' / ^ 

Finally, from Eq. (|25|), G can be written as 2s^(l -|- s^H{s,a)^ with H{s,a) uniformly bounded by 3. Then, 

V2G(s, a) = 2s, (1 + slH{s, a)) (28) 

Gathering Eqs. (p6|),(27),(p8|), we get the real-analyticity of ^^'^"-^ for {s,a) G R x /, where / = {—ao,ao)i for 
some aQ > small enough. 

■ 

We next compute /o, /i, /2 and obtain the following proposition. 

Proposition 3 For a small enough, the function , e = ± defined above and its derivative with respect to s 
have the following expansion 

X'{a,r,s) = /J(s,r) + /f(s,r)a + /|(s,r)a2+0(a3), (29) 

§-,X%c.,r,s) = ^J^(s,r) + -^Jf{s,r)a + ^Jiis,r)a^+Oia^) (30) 

where ff , ^ = 0, 1, 2 are defined as in ^ and \0{a^)\ < C\ a I , with the constant C independent o/ s G M and 
rG [0,1). 
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Proof of Proposition ^. Wc will prove the proposition only for ■ To proceed, we list the expansions of 
the form (|2^) for several quantities, obtained after elementary computations. 

i^{a,r,s) = TT - s - TTs^a + 2(1 - r)ssa^ + O(a^) (31) 

2 

6'(a, r,s) = 7r-2a(l+r) + ^a2^0(a3) (32) 

—oiSs 1 — ^Q;^ 

Z-{v{s)) ^ \ aSs 1 +0{a'). (33) 



-l + ^a^ 



Using Eqs. (|l|),(|3|), we get that 



3in (V'(a,r,s)) = + 7r.s,c,a - (2(1 - r)s,,c^ + ^s^ja^ + O(a^) (34) 

2 



COS (V'(a,r,s)) = - Trs^a + (2(1 - rjs^ - ys^c^)^^ + C'(a^) (35) 

2 

sin (6'(a,r,s)) = 2a(l + r) - + ©(a^) (36) 



cos 



{e{a,r,s)) = -l + 2a2(l + r)2+C)(a3). (37) 



Using Eqs. (DJ), (34) and (M) we obtain 



where 



(—Cs + TTSgOf — Sg — TrCsS^a —(1 + Cs)a\ 
Ss+^SsCsa Cs+TTsla s^a \ + n{s)a^ + 0{a^) , (38) 

-(1 + Cs)a -s^a 1 / 

/c, + 4c.s' + 1-2(1 + r)s2 ^ + 2c,s,(l + r) + 4s3 tts^ 
^('S)= -f-2c,S3(l+r)-4s' -2(l + r)s2 + 2^c,s2 ^s^c. 



and using Eq. (p^), we have 



'^ct^ct(l — ^s) 

e'^+N= I -s„s, I = I -ass I +0(a^). (39) 




An easy computation yields 





Zl{vis))={ -c2(s) 6(s)c(s) I = ( -a^s^ as, |+0(a^). (40) 
fe(s)c(s) -fo2(s) 




Using Eqs. (ha), (36) (BTp and (p9|) and the previous equation, we get 



e 



a(l + 2r + c,)-a^f 

e(a,r,.)z_(.(.))gsx+^^ [ I +0(a3). (41) 

-1 + (l + 2 r + 2 + Cs + 2 r Cs) 



Applying e''^^*'^- to the previous equation and using Eq. (p8[), we finally get Eq. ( p9| ) for e = +. The expression 
of the derivative ( pO| ) is then an immediate consequence of the analiticity of the function x'^ ■ 



The results for x a-nd -^x a-i'c obtained similarly together with Eq. (|1^ 



Since the quantities of the form 0{a^) in Proposition || satisfy |C'(q;'^)| < Cja"^! for some C independent 
of r, the expressions (||)-(||) are straightforward consequences. Hence the shape of the extremal front at time 
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T = /ca/tt is known for a small. In particular its image with respect to the map defined in Section |i| is 
approximated, in the sense, by a circle of radius 2r(a) centered at the origin. 

We finally note that, for fcjv/ even. Lemma ^ is still valid, while, with computations similar to those made 
in the proof of Proposition ||, it is easy to see that the formulas for /|, fc = 0, 1, 2 simply differ, with respect 
to (^, for the sign of the first two components. 

Remark 6 Repeating the previous computations, we also obtain series expansions for £'^{s,kM — and 
■^£'^{s,kM — l,a). Indeed, we just have to replace r hy 1 + r. In that case the shape of the extremal front 
J^{a, {kM — 1)7''), after applying the map M„, is approximated, in the sense, by a circle of radius 2(1 + r(a)) 
centered at the origin. 



4 Case r{a) = r G (0, 1) 

In this section, we study the case in which a tends to zero with r{a) — f, for a constant f G (0, 1). More 
precisely we consider the decreasing sequence ~ 2(k+f) ' ^'^'^ k> \. We first describe the minimum time front 
at T = fcjv/Tr, then we identify and study the candidate for the limit synthesis and finally we prove Theorem |l|. 

4.1 Description of the minimum time front at T = /c^/tt 
The purpose of the paragraph is to prove the following proposition. 

Proposition 4 Fix 5 > Q small. For a small enough with r{a) > 5, the extremal front !F{a^ /ca/tt) is home- 
omorphic to a circle. As a consequence, the switching curves defined inductively in Eq. are optimal up to 
k = kM and OF(a, /cm""), the minimum time front at time /cm"" coincides with J-'{a, ki\iTT). 

Proof of Proposition ^ From Proposition ^ we get that the extremal front J- {a, k^Tr) is the union of two 
arcs, £^{a,s), s G [0,7r] and £~{a,s), s € [0,7r] so that, for e = ± and s G [0, tt], 

'2r{a)eaCs 

f^(a,s)==l 2r{a)eas, l+Oia"^), (42) 



and 



-1 



s) = 2r(a)ea | c, | + 0{a^). (43) 




Moreover, at s = and s = tt, the derivatives of £^(a, s) are only one-sided, i.e. as s > tends to zero and s < tt 
tends to tt. By a trivial continuity argument, one can parameterize J^{a, fci/Tr) as a closed continuous curve 7 
defined on [0, 2n] so that 7(5) = £~^{a, s) for s G (0, tt] and 7(5) = £~{a, s — tt) for s G (tt, 27r]. Moreover, with 
the previous computations, it is immediate that 7 is in fact piecewise with possible discontinuity jumps for 
^7 at s = and s ~ tt. 

Since the curve 7 is in a neighborhood of the south pole of size proportional to a (thanks to Eq. (p^), it is 
enough to prove that the orthogonal projection 71 of 7 on the [xi, a:;2)-plane is homeomorphic to the circle e'^" , 
s G [0, 2tt]. Using Eq. (^), we see that ||7i(s)|| = 2r(a)a + 0(p?) on [0, 27r], which implies that the continuous 
function ||7i(s)|| is always strictly positive for a small enough. We can therefore parameterize 71 using polar 
coordinates {p,P), i.e., for s G [0, 27r], 

j^is) ^ pis)e^^^^\ 

where p(-) :~ ||7i(-)ll ^-^d the function /?(•) are defined on [0, 27r], continuous and piecewise C^, with possible 
jumps of discontinuity for their derivatives at s = and s — n. 

In addition p(0) = p{2n), /3(0) ee (3{2n) = tt (mod 2tt) and, from Eq. (||), P{s) = tt - s + 0{a). To prove 
Proposition^, it suffices now to prove that /3 is a monotone bijection from [0, 2n] to [— tt, tt]. The latter simply 
results from Eq. (^). Indeed, from that equation, we get that ^/3(s) = —1 + 0(0;) where /? is differentiable and 
the one-sided derivatives at s = and s = tt verify the same equation. We deduce that /3 is strictly decreasing 
for a small enough. 
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Wc next show that OF{a,kM'^), the mimmum time front at tmie fcA/7r coincides with J-{a,kMT^)- By the 
resuhs of |^], we first notice that any time minimal trajectory starting at the north pole reaches the south 
pole in time T > kMT^- Therefore OF{a, /cmtt) is not empty and is included in ^(a, kMi^) according to the 
PMP. According to Theorem 27 of OF{a,kMT^) is a one-dimensional piecewise compact embedded 
submanifold of S^. By an easy topological argument, we deduce from the above that OF{a,knjTT) coincides 
with J^{a, fcjvfT'')- B 

Remark 7 Thanks to Remark ^, and with arguments similar to those of the previous proof, one can prove that 
J^{a, — l)7r) is optimal for a small enough, with no assumptions on the remainder r. 



4.2 Optimal synthesis for the hnear pendulum control problem 

Recall that Ma : R"^ — > is the composition of the projection {xi,X2,X3) i-^ {xi,X2) followed by the dilation 
by 1/a. With the results of the previous subsection, it is clear that the original control problem on can be 
reduced, near the south pole, to a planar control problem on the neighborhood of the south pole delimited by 
OF{a, km"^) ■= Ma{OF{a, knj'n')) along (S)q, the control system obtained as the image of by Ma, i.e. 

(y^ . r = -cos(a)z2, , . 

^ 1 Z2 = cos(«)zi + - (azi)2 - (az2)2, (zl,Z2)eR^ \u\<l. ^ ^ 

It is therefore natural to conjecture (simply set a = in OF{a, /cmtt) and (S)q) that the Hmit synthesis should 
be that of connecting the circle of radius 2r(a), C{0,2r{a)), to every point of the disk 5(0, 2r(a)) along the 
control system {Pen) given by Eq. (^, which we rewrite as 

{Pen) z = Aoz + ubo, with Aq = , ^" ^ (l) ' ^'^^^ 

where z G R^ and u e [—1,1]. The control system {Pen) corresponds to a linear pendulum with a forcing term. 

Theorem ^ simply states that the conjecture is correct and, as a first step for an argument, we describe, in 
more details in this subsection, the conjectured limit synthesis. Hence we focus on the following problem. 

(P) Fixed p s]0, 2], for any given y £ -6(0, p) Gnd a time optimal trajectory connecting the circle of radius 
p centered at the origin to y along the control system {Pen) . 

Remark 8 The problem of computing the TOS for the linear pendulum, taking the origin as a source, is studied 
in any textbook of optimal control. Here as source we take the circle of radius p centered at the origin, which 
is a level set of the Hamiltonian H = \{z1 + z|) associated to the uncontrolled system. 

It is easy to see that the solutions of problem (P) must be bang-bang trajectories. Indeed since {Pen) 
is a bidimensional linear control system it is well known that this property is guaranteed by the Kalman 
controllability condition det(6o,Ao6o) 7^ 0, which is satisfied by {Pen). To determine the TOS, we first look 
for the switching curves. We know that every extremal trajectory for the problem (P) must satisfy the 
transversality condition of the PMP stated in Section Here the source manifold is the circle C(0,/5), and 
the transversality condition essentially translates into the property that the vector A(0) (that, without loss of 
gencrahty, we will assume unitary) is proportional to z{Q) S C(0,p) (identifying the cotangent space with the 
plane R^). To determine completely A(0), it is enough to observe that a necessary condition for z{.) to be 
optimal is that i(0) points inside the disk i?(0, p), i.e., if we denote by Uopt the corresponding control, then 

< z(0), z(0) >< ^ < 2(0), ^o^(O) + uoptho >< < z(0), Uopt6o >< . 

Therefore, Uopt = — sgn< z(0),6o >• On the other hand, from the maximality condition of the PMP, we have 
Uopt =sgn< A(0),6o > and, therefore, one can define A(0) —z{0)/p. Finally Uopt = ^sgn(z2(0)) (except at 
the points ±(p, 0) ), while the switching time tsw must satisfy the condition < A(tstu), ba ^■2{tsw) = 0. 
Consider now the adjoint system 

i : (46) 

A2 — Ai. 
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If we identify R"^ with the complex plane, so that z ~ zi + iz2 and A = Ai + iA2, then the equations (|45|), ( p6| ) 
become 

i = i(z + u) and A = lA. 
Moreover we can set z(0) = — pe^'^ and A(0) = e^'^ for some 9 E [0, 2tt[ and the corresponding solutions are: 

z{t) = (z(0) + Uopt)e'' - Uopt = -pe^^*-") + Uopt{e'' - 1), 
X{t) = A(0)e'* = e'(*-''). 

The switching curves are determined by the relation tg^ = 9 (mod tt) and this allows to conclude that the 
switching curves are the following two semicircles of radius 1: 

z{9) = 1 - p - e''^ 6'e[0,7r[, 
z(6i) = p-l-e''' 9 e [tt, 2tt[. 

These switching curves cannot be opti mal for p < 2 since they are not locally optimal, as can be easily checked 
using the definition given in Section 2.2, We conclude that the optimal trajectories are bang arcs and the 
corresponding control depends on the sign of the component 2:2(0) of the starting point. 

To conclude the description of the synthesis, it is enough to determine the cut locus, i.e. the set of points 
that are reached by two or more optimal trajectories at the same time. Assume that 2 € C belongs to the cut 
locus. Then, there exist s € [0, tt), s' G [tt, 2tt) and t such that 

z = -pe*(*-^') - 1 + e". ^ ^ 

Therefore |z — 1 + e**| = |z + 1 — e**| = p. In particular, denoting by z the complex conjugate to z, we have 

(z - 1 + e**)(z - 1 + e~**) - (z + 1 - e'*)(z + 1 - e"**) = -4zi + 4zi cost + 4z2 sint = 0, (48) 
(z- 1 + e**)(z - 1 + e"'*) + (z + 1 - e'*)(z + 1 - e"**) = 2zf + 2z| + 4 - 4cosi = 2p^ (49) 

From ( ^) we have that cost = , and, substituting in (^), we find that z must satisfy the equation 

zt + z^ + 2zlzl - p^zl + (4 - p^)zl = 0. (50) 

The previous computation show that the cut locus is a subset of the set of points belonging to the locus defined 
by (^). Actually it is easy to see that this is the proper subset obtained with the additional condition Z1Z2 > 0, 
that corresponds to t < tt. The precise shape of the optimal synthesis, which is now clear, is portrayed in Fig. ^ 
for a particular value of p < 2. Notice that, from the previous computations, wc have pe** = pe*^ + 2 — 2e** 
and, since pe'^ + pe*'* = 2pe'* — 2 + 2e** and pe'" — pe*'* = 2 — 2e'* are orthogonal in the complex plane, we 
find easily the following equation: 

(2 — pcoss')(cosi — 1) — psins' sint = 0. 

Consequently, for t G [0, 27r[ and s' G [tt, 27r[, one has, along the overlap curve 

. ,-. psins' 

t = t(s') = -2arctan— ^- . (51) 

^ ' 2-pcoss' ^ ' 

This expression will be useful in the following. Also, notice that combining ( ^ ) and (^l|) one easily finds a 
parametrization of the overlap curve in terms of s' and that in an analogous way it is possible to parameterize 
it by means of s. From now on we will denote by 7pen(') the parameterization of the overlap curve with respect 
to the parameter s. 

Remark 9 If p = 2 the previous reasoning does not apply and indeed the synthesis is different. In this case the 
overlap curve coincides with the switching curves and with the trajectories reaching the origin corresponding 
to u = ±1. A simple way to prove this fact is to study the optimal synthesis starting from the origin with 
vector fields with opposite signs, and observe that the extremal front at time tt is a circle of radius 2. 
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Figure 7: The overlap curve for the pcnduhim problem 




Figure 8: Propositions || and || 



4.3 Proof of Theorem |l] 

The proof of Theorem |l| is divided in two parts. Roughly speaking, defining := Ma[C\^^(^)) for e = ±, 
we will look separately at the shape of the synthesis far from P" and inside neighborhoods of P", e = ±. Let 
us call Vq, the image with respect to Ma of the neighborhood of the south pole enclosed by OP(a, /ca/tt) and 
the ball of center P" and radius ^. 
Then, the previous cases correspond to the following two propositions whose meaning is clarified by Figure ^. 

Proposition 5 Let f G (0,1) and a^. :~ 2(k+f) • Then, for any ^ > there exist a positive integer k and a 
compact interval I C (0,7r) such that it is possible to find a curve defined on I for k > k, verifying the 
following: 7^ divides V^^ \ (P"'°(^) U P"'°(C)) in two connected components, such that "above" 7^ the optimal 
feedback associated to the synthesis for a = oif^ takes the value —1 and, below 7J;!, it is equal to 1, and in 
particular 7^ is an overlap curve for a = a^. Moreover, 7^ converges to Jp^n ^"^ topology of I. 

Proposition 6 Consider the notations defined above. Then there exist ^ > 0, r^, e = ± with < r_ < t+ < tt 
and a positive integer k such that, for every k > k, it is possible to find two curves 7° j, and 7^J_ j., defined 
respectively on [0,r_] and [T+,7r], verifying the following: 7°^, divides Va^f^B"'' {^) in two connected components, 
such that "above" 7° j, the optimal feedback associated to the synthesis for a — takes the value —1 and, below 
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7° J, , it is equal to 1 , and in particular the 7° j, are overlap curves for a — ak- Moreover, 7° ^ and 7^ converge 
to 

Ipen *^ ^'^s C'*' topology respectively o/[0, r„] and [T+,7r]. 

The choice of studying the synthesis separately in neighborhoods of and far from is justified by the 
fact that the proofs of the previous propositions rely on different implicit function arguments. 

It is clear that, combining Proposition ^. for an appropriate choice of ^, with Proposition ^, one almost 
completes the proof of Theorem |l|. We will not prove explicitly that the convergence of 7^ to 7°g„ with respect 
to the parameter f is uniform in any closed interval of (0, 1). As explained in Remark nQ. this can be done 
with the same methods used in the proofs of Propositions |^,^, but with much more computational efforts. 

We will therefore only provide the complete proofs of the propositions. For this purpose, the first step 
consists in checking whether the switching curves Cf. , e = ±1 are optimal or not. In that regard and 
similarly to the case of the linear pendulum, we have the following result: 

Lemma 3 Let r g]0, 1[. Then, if a is small enough and r{a) = f, the switching curve Cf. is nowhere locally 
optimal i.e. all the extremal trajectories switching on Cf.^^ lose optimality before reaching it. 

Proof of Lemma |^. For simplicity we define S{s) := C'fejj^l-?) a-i^d we assume kM odd. As in the proof of 
Lemma ^, we get the following asymptotic expansions, after applying the map A/q,: 

Sis) ^ ['^-l;^^)+Oial .(0).(2^7(")), (52) 

S'is) = [-:;)+Oia), + (53) 

S"is) ^ )+0(a). (54) 

Integrating the above equation, we have 

S'W = S'(0) + y^'s"(r)*=(-;'+^(^;>'), (65) 
.(.) . 5,0) + /-S'M...(^-J;- + °W) (66) 

and therefore 

J_X^(S(,))=f ~^2(^) ^Ss + Oisa) 

^^^±{^{S)) ^ s,is)±'-^^l-a^S,is)^-a^S2is)^ J { 2r ~ 1 + c, + Oia) ± il + Oia^)) 

Here Si, i = 1,2, denotes the i-th component of S. Dividing the above equation by 1 + Oia), we can assume 
that the first component is identically equal to —Sg. The same can be done with the expression (p5|), so that 
it is possible to compare the three vectors obtained in this way simply by looking at the second components, 
which are equal respectively to 2r — 1 + Cs ± 1 + C'(a) and Cg + Oia). In particular, the fact that £'(•) is nowhere 
locally optimal if a is small enough follows from the inequalities 2r — 2 + Cs + Oia) < Cg + Oia) < 2r + Cs + Oia). 

m 

A straightforward consequence of the previous result is the presence of a non trivial cut locus in the 
neighborhood of the south pole enclosed by P(a,fcA/7r). It remains to clearly define that cut-locus, which is 
the purpose of Propositions || and ^. 

4.3.1 Proof of Proposition 1^ 

As usual, we only provide an argument in the case km odd and we fix the remainder equal to f £ (0, 1). 



Recall that, according to Section 4.1, OFia,kM'^) is approximately (up to order a^) a circle of radius 
2m. To describe the synthesis inside the neighborhood of the south pole enclosed by OF(a, /ca/ti"), it is more 
convenient to use the two dimensional control system which is rewritten as follows by using Eq. (H), 



z = CqAqz + u— -v/l — a^llzP&Q. 
a 
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Wc set X^{z) := CaAoz + e^y/l - a'^Wz^bQ and ^^^(z) := Aoz + ebo, for e = ±, and wc define OF{a,kMi^) 

as the image by Ma of OF{a, kMT^)- Then we know that, up to order a, OF{a, kMT^) is a cirele of radius 
2r. In particular, as in the proof of Proposition one ean eonstruct a piccewise smooth parameterization 
ctq : [0, 2tt] F{a, /cm tt) so that 0-0(0) = P", ^^(Tr) = P" with a loss of regularity only occurring at s = 0, tt 
(with two-sided differentials at any order). In particular (Tq,(-) approximates in the sense the function 
a : [0, 27r] ^ C R^, defined as cr(s) = 2fe^^^~^\ which is a parameterization of the circle of radius 2f. 

Taking into account Lemma |[ the cut-locus in Va is contained inside the set of points Q G , besides , 
such that there exists (s,s',t) € (0, tt) x (tt, 27r) x (0, 27r) for which Q = e*'^+CTQ,(s') = e*^-aa{s). 

In view of applying an inverse function result for characterizing this set, we consider the map $ defined on 
[0, tt] X [tt, 2tt] X [0, tt] by 

$(s,s',i) := (s,e*^r"a(s')-e*^'-''V(s)), 

which takes values in R"^. Similarly, for k > 1 and Uk as in the proposition, we consider the map defined on 
[0, tt] X [tt, 2tt] X [0, tt] by 

$fc(s,s',t) := (s,e*^+V„,(s')-e*^-'^a.(s)). 

Note that, since the vector fields X^'' converge uniformly to X?"^" on Va, it is easy to see that $fc converges to 
$ in the norm. 

For (Pen), a point of the overlap curve, besides P", is then identified with a triple {s,s',t) S (0,7r) x 
(tt, 27r) X (0, tt) such that $(s,s',t) = (s,0,0). In other words, the overlap curve can be parameterized by 
means of the map w : [0, tt] R"^ defined imphcitly by $(w(s)) = (s,0, 0), while 7°g„ can be obtained as the 

composition of the two maps e^'^- <t{s) and w{s). 

Similarly, we would like to define the overlap curve corresponding to (S)^*. , for k large enough, by means of 
the function Wk defined by ^kiwkis)) = (s,0, 0). To proceed, we will apply Theorem ||. The first task consists 
of computing det along the overlap curve. 

J -y-pen , -O-pen 

Lemma 4 Along the set of triples (s,s',t) G (0, tt) x (tt, 27r) x (0,7r) for which e + cr(s') = e - o-{s), we 
have 

4r f 1 — ) sin s' 

detV^{s,s',t) = ^-i — ' ,_ . — —T. 

(1 — r cos s')^ -I- (r sm s'Y 

Proof of Lemma ^. One has 

da 

' ds 



det ^$(5, s^ t) = det ((e*^r" ), If "e**r"a(s') - l^^"e*^"V(s)) . 



By taking into account that $(s, s',t) = 0, the previous determinant is equal to twice the first component of 
(e*^?'"),|^, i.e., detX"I>(s,s',t) = 4fsin(s' - t). Using Eq. (js^), one concludes. ■ 

Observe that detD$ ^ if s' ^ 0, tt. In particular, if we consider a closed interval / C (0, tt), then the 
set Im(w)|/ plays the role of the compact set K in Theorem ^. All the assumptions of the theorem are then 
verified, and therefore we have proved the existence of a map Wk defined on /, satisfying ^k{wk{s)) = (s,0, 0) 
and converging uniformly to w. If we define 7^ as the composition of the two maps e^'^- cr(s) and w{s), then, 
since I was chosen arbitrarily, the proof of the theorem is complete. 



4.3.2 Proof of Proposition § 

With the previous notations, let ipk be the smooth map defined on [0,7r] x [tt, 27r] x [0, 27r] by 

(Pfc(s,s',t) = e'^+'ak{s') - e'^-" (Jk{s). 

For the rest of this paragraph, we drop the index k to get lighter notations. 

From the Taylor expansion of Lp around the points (0, 27r, 0) and (tt, tt, 0), we derive the asymptotic behaviors 
of the cut locus close to the points P", e = ±, since that cut locus belongs to the level set = 0. We will only 
perform computations at (0, 27r, 0) since they are entirely similar at (tt, tt, 0). 
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Let us call the two components of ^p. We use ipi' to denote the partial derivative of the component 

(^'*^ with respect to s evaluated in (0, 27r, 0) and we define in an analogous way all the (multiple) partial 
derivatives evaluated in (0, 27r, 0). Set s := s' — 27r. Then, after computations, we have ip^p = ip'p = ip"^^ = 
and 

^i^; ^ ^2f + 0{a) ^14' = 2f + 0{a) , - 2 + 0{a) , 



(^i^' = 2f + C'(a), (^f ' = -2f + 0(a) , (^f ^ = 2f + ©(a) . 



We thus get 



= -2fs^ + 2rS^ + 2f -4:fst + 4:rSt + 0{a\{s,S,t)\'^) + 0{\{s,S,t)\^), (57) 

and 

</7(2)(s, t) = ,^(2)5 + + ^(2)^ + 0(1(5^ i)|2) ^ 2fs - 2rS ^ 2t + 0(a|(s, I, t)|) + 0(|(s, I, t)|2) , (58) 

where, here, ©(•) is imiform with respect to a. 

Fix ^0 > small. We are looking at the cut locus in a neighborhood of P", and thus, we can assume 
|(s, s,t)\ < ^0 for some ^0 > 0- The purpose of subsequent computations consists of expressing s < and t > 
as functions of s, for < s < ^0, by using the equations </j^^^ — and (/j'^^) = q. 

From (^(2) = 0, by applying the implicit function theorem, for ^0 small enough and |(s,s)| < ^0, we get 
t = h{s, s) with h G C^. Moreover, since h{s, s) = 0(1(3, s)\) we have 

h{s, S) = rs-rs + 0{a\{s, §)\) + 0{\{s, s)\^). (59) 

Consider now the map 

S) = ; , 

s — s 

which is well defined and for s > 0, s < 0. Again, it is possible to apply the implicit function theorem to 
the equation = 0, so that we get s as a function of s and this gives the existence of the overlap curve. 
Moreover, by combining Eqs. (pq) and (p9h, we get the following 



0(s, §) = ,s(l + f) + S(l - f) + 0(a|(s, + 0(|(s, I)|2) = 0, (60) 
and then, \s\ = 0{\s\). Therefore, from this estimate and the above ones, we immediately obtain that 

1= -(ii^ + 0(a))s + 0(s2), (J!_+0(a))s + 0(s2), (61) 

from which we get that the overlap curve converges uniformly to 7pen- 

The proof of Proposition |^ is now complete. I 

Remark 10 In order to prove that the overlap curve 7^ converges to 7pe„ uniformly with respect to f in any 
closed interval / C (0,1), it is enough to follow the lines of the proofs of Propositions ||, ^ by considering 
f as an additional variable. For instance, for the Proposition ^, one needs to define the maps $ : — * R'', 
$(f, s, s',t) := (f, s,e*^r"tT(s')-e*^-"V(s)) (recall that cr(-) depends on f) and $fc : ^ K^*, $fc(f, s,s',i) := 
(f, 3,6*^^+ CTfc(s') — e*"^- CTfc(s)) (where at = 7r/(2(f + fc)) and (Tfe(-) depends on f). 

The uniformity with respect to f is then proved by applying Theorem |^ to $, $fc with /C = {$~^(f , s, 0, 0) : 
(f, s) G I X J} where / x J is a compact subset of (0, 1) x (0, tt). 
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5 Case r = Ca 



5.1 Description of the minimum time front at time fcj\/7r 

Fix C > and consider the sequence (a.k) such that r(afc) = Coik, k > 0. As before, we drop the index k when 
possible. For ak small enough, one deduces, from the analysis of 0, that the south pole is not reached at time 
fcmTT = [^^l""! so that the optimal front at time fcM'"' is not empty. The next result provides a description of 
the extremal front at time fcj\f tt. 

Lemma 5 Define the planar curve C : [0, 2tt] by 

,s«(7r + 2C-^.s2/2) 
Then, for s G [0, n], we have 

£+{a,s) = {a'^C{s),-lf + Oia^), (63) 

and 

^£+{a, s) = {a^^£{s), Of + 0{a^). (64) 
OS as 

At s = and s = tt, the derivatives are only one-sided, i.e. as s > tends to zero and s < tt tends to tt. 
Similarly, we have, for s G [0,tt], 

£-{a,kMTT,s) = {a^£{s + TT),-lf + 0{a^), (65) 

and 

kMTT, s) = + ^), 0)^ + 0{a^), (66) 

with one-sided derivatives at s = and 3 = 11. 

Proof of Lemma |5|. This is immediate from Proposition^ applied in the case r(a) = Ca. ■ 
For C < 7r/4, consider 9d G (0, 7r/2) with sin{9d) = jn. The curve £(s) has two double points 
D"*" = £(5]*^) = £(sj), with s\ = Qd and S2 = tt — 9d, and D~ ~ £{si ) — £{82 ), with s^ — tt -{- 9d and 
S2 ~ 2tt — 9d. It also has four cuspidal points Cpf , i = 1,2 and s — ±, corresponding to the values s = s^uisp 
where sin^ s = 2±^£/£_ 

Finally, let a be the closed Jordan curve defined as the restriction of C{s) to [0, s^^] U [sj, sj"] U [s2 , 2tt]. If 
C > 7r/4, we simply define a to be £. 

At the light of the previous result, we get that J-{a,kMTT), the extremal front at time fc^/Tr, is contained 
inside a neighborhood Wa of the south pole of order 0{a^) neighborhood of the south pole. Therefore, in order 
to understand the shape of the optimal synthesis inside Wa , we must rescale the whole problem by Na , the linear 
mapping from to R'^ onto defined as the composition of the orthogonal projection {xi,X2tX^) s- {xi,X2) 
followed by the dilation by l/a^. 

For X G Wa, we first consider (A)q,, the image of (E) by Na, i.e. (A)^ is the planar control system given by 

^"^^^ '■ { 1 = CaZ^ + u^^l-a^z^. ^^^^ 

Let Ca be the image of T{a, fcj\/7r) by Na- From Lemma ^, Ca converges to C in the topology. It is clear 
that, for C > tt/4, Ca ■ [0, 2tt] is homeomorphic to e'**, s G [0, 2tt]. In the case where C < tt/4, the next 

lemma shows that, for a small enough, Ca has the same shape as C. 



Lemma 6 If C < tt/4, then Ca is described by the following picture, where Cpf{a) = Ca{s^cusp i — ^,2 

and e = ±, are cuspidal points and D'{a) are double points with 

D+{a)^Ca{st{a)) = Ca{st{a)), {a) = Ca{s^ (a)) = Ca{s2 {a)), (68) 

where s^^^p j(a) and sf{a) tend respectively to s^^^p ,^ and sf as a tends to zero, for i — 1,2 and e = zt. For 
a small enough, set Oa, the closed curve defined as the restriction of Ca{s) to [0,sf{a)] U [s2 (a) , s^ (a)] U 
[s^(a),27r]. Then, it is a Jordan curve. 
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Proof of Lemma |6|. For i = 1,2 and e = ±, the existence of the cuspidal points Cpf{a) is obtained by 
applying the implicit function theorem to the equation DL{s,a) = 0, where the function DL{s,a) is defined 

by ^ 

DL{s,a) := — 
as 

in the neighborhood of each (s^w^p ii 0)- We have 

and we conclude. The uniqueness of these four points, on [0, 2tt], is trivial since DL{s, a) = j^C{s) + 0{a). 

Similarly, for e = ±, the existence of the double points D^{a) follows after applying the implicit function 
theorem to the equation DP{s, s' , a) = 0, where the function DP{s, s' , a) is defined by 

in the neighborhood of each (sf , Sj, 0). For the uniqueness, we proceed as before. ■ 
In the case C > 7r/4, we also define ctq to be equal to Ca- As a consequence, we are able to characterize 
OF{a, kMT^), the minimum time front at time /ca/tt when C =^ tt/A. 

Proposition 7 For a small enough and C ^ 7r/4, the minimum time front at time k^TT, 0_F(q;/£, fcjv/7r) is 
equal to CFa, the inverse image on S'^ , by Na, of Ga. 

Remark 11 As a consequence, we deduce that, for C > tt/A and a small enough, the optimal synthesis between 
J-{a, {kM — l)7r) and J-{a, fcj\/7r) is simply given by the extremal flow whereas, for C < 7r/4, there is a loss 
of optimality along certain extremal curves starting at J-{a,{kM — l)7r) before reaching J-{a,kMT^)- The 
values of s corresponding to such curves can be deduced from the previous characterizations of J-'{a, kj^i-n) and 
OF{ak,kMTT). 

Proof of Proposition ^ Recall that OF{ak, fcA/Tr) is a piecewise submanifold of ^{a, kMT^)- As in the 
proof of Proposition ^, the result to establish is a consequence of the fact that ~ Ca '■ [0, 27r] M-^ is 
homeomorphic to e*'*, s G [0, 27r] and it can be achieved by means of simple topological arguments. 

In the case C < 7r/4. a a is a piecewise Jordan curve homeomorphic to e*^, s € [0, 27r]. A simple 
topological argument yields the conclusion. ■ 
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5.2 Limit of the synthesis 



It remains to describe the limiting dynamics close to the south pole. In order to take the limit, as a tends to 
zero, in (A)q,, one must reparameterize by the time at. The limit is then given by the control system 

(A) : ( = °' 

[ Z2 = U. 

We now describe the optimal synthesis for the limit problem, i.e. for the problem of reaching in minimum time 
every point inside a along (A) and starting from a. Because of the symmetries of a and because the tangent 
vector to a is vertical only at s = and s = tt, there exists a unique overlap curve {Seg)c, defined as the 
segment of the ^i-axis between the points (— 2C, 0) and (2C, 0). Above it, the input u takes the constant value 
— 1 and, below that overlap curve, the constant value 1. Integral curves are clearly vertical lines. 

We next intend to prove that the optimal synthesis consisting of reaching in minimum time every point 
inside ctq along (A)^ and starting from ctq, converges to the previous synthesis in the following sense. 

Theorem 2 Assume that C ^ i: j/^. As a tends to zero, the time optimal synthesis associated to (A)q, inside 
Ga is characterized by an overlap curve {Seg)'^, converging to {Seg)c in the topology, and, above {Seg)^, 
the control u takes the constant value —1 and below [Seg)^, it is equal to 1. Moreover, there exist only two 
time optimal trajectories reaching the origin and, in the case C Kir/A, these trajectories start from D^, £ — ±, 
the double points of Ca ■ 

Proof of Theorem |2|. Fix C ^ 7r/4. We first notice that, for a small enough, there are not switching 
curves inside CTq. Therefore, the cut-locus may only occur as images by iV„ of points M G S*^ such that 
M = ei'^-5-{s) — ei^+a{s') for t e [0,^]., s e [0,7r] and s' e [tt, 27r]. Proceeding exactly as in the proof of 
Theorem^, we apply inverse function arguments first in neighborhoods of (7q,(0) and (7Q,(7r) and second in a 
region enclosed by ctq excluding such neighborhoods. It is then easy to determine the values of the input u in 
each connected component of the region enclosed by ctq minus (Seg)'^. 

By a continuity argument, it is clear that there exist only two time optimal trajectories reaching the origin, 
one above {Seg)'^ and one below. Finally, suppose that C < n/A. In that case, it was proved in that the only 
extremals starting at a point £q(s) and reaching the origin from above the overlap curve {Seg)^ correspond to 
values of s verifying one of the following three possibilities as a tends to zero: (a) s tends to zero, (b) s tends 
to 7r/2, (c) Ca{s) is a double point also associated to s' = v{s) — s. In view of what precedes, only possibility 
(c) is allowed for optimality. Theorem ^ is proved. ■ 

Remark 12 As a consequence of the previous argument and from the results of we get that, for a small 
enough and C < 7r/4, 

s^(a) = v{s^{a)) — s^{a), S2 {a) = 2tt + v{s^ (a) — tt) — s^ {a), 
where sf (a), i = 1, 2 e = ±, were defined in (|68|). 



6 Case r(a) = 



We assume here that r{a) = 0, i.e. ak = ^jr for fc > 1. From Proposition |^, we know that the extremal front at 
time ([^] — l)7r = 5^ ~ ""j encloses the south pole, is optimal and is approximately (in the sense) a circle of 
radius 2r{a)a around the south pole. Moreover, at time [^]7r, we know that the extremal front must contain 



the south pole and is equal, up to O(a^), to {a^C, —1)^ given in (63) and ( |65| ) with C = 0. In that case, the 
minimum time front reduces to the south pole. 

In this case it is interesting to consider the synthesis starting from the extremal front at time (fc^/ — and 



it is natural to compare it with the synthesis of the linear pendulum studied in Section 4^ and corresponding 
to p = 2. Let us first describe briefly that synthesis. Let D2 and C2 be the disc and the circle centered at the 
origin and of radius 2 respectively. The overlap curve inside D2 coincides with the switching curves and with 
the trajectories, corresponding to u ~ ±1, connecting the points (±2,0) to the origin. In particular, it means 



25 



that an optimal trajectory of the synthesis starting at any point P G C2 reaches the origin, and thus, there 
exists an infinite number of optimal trajectories from C2 to the origin. 

For a > and r{a) = 0, the situation is rather different. Let us first define F{a, {kM — f )7r) to be the image 
of JF(a, {kM — l)7r) by Ma- Then, for a small enough, it was shown in j^, that the only optimal trajectories 
starting from F{a, {kM — l)7r) and reaching the origin are those starting at P" and P". Let us refer to them 
as 7+ and 7". Therefore, in the case r{a) = 0, the synthesis for a > is rather different than the synthesis of 
the limit candidate when a tends to zero. It is a clear indication that the case r{a) = is more delicate than 
the cases r{a) positive constant or r{a) = Ca. However we are still able to give a partial description of the 
limit synthesis as the next proposition shows. 

Proposition 8 Assume that r{a) — and a is small enough. Then the switching curve (resp. ) 

is optimal for some interval [0, s(a)], s{a) < n, and it is above (resp. below) 7+ (resp. j~ ) as long as it is 
optimal. Moreover, we have 

lim s{a) = s := arccos ■\/l/3. (69) 

a — ^0,r{Q)— 

Proof of Proposition |^. We only provide an argument for , being the other case analogous. To prove 
the first statement of the proposition, we reason by contradiction. If the switching curve is not optimal on any 
interval [0,t], t > 0, we get the existence of an optimal trajectory starting at !F{a, {kM — ^)'^) above P+ and 
reaching the origin, which is equal to the concatenation of an integral curve of X- and a piece of 7^. Therefore, 
an optimal integral curve of starting above 7"*", must either switch or lose optimality before reaching 7+. 
If the second possibility occurs, we must have an overlap i.e., at that point an optimal integral curve of X+ 
arrives. Close to P+, the latter would imply that the optimal integral curve of X+ starts at ^{a, {kM — l)7r) 
above P^. This is impossible because, from every point of £~^{a, {kM — ^)'^), the value of the optimal control 
is —1. Let s{a) < tt be the first value of s for which ceases to be optimal. Define 

dC^ 

H{s) det {X+{C+Js)), ^(s), (s)) , 

for s G [0, tt]. Then, s{a) is the smallest solution in (0, tt] of H{s) = 0. It is easy to see that H must take the 
value zero before tt. We deduce that s{a) < it. By taking the asymptotic expansion of the previous expression 
as a tends to zero, we get 

H{s) = '!^s,a^{l + 2,cos{2s) + 0{a)). 
Then s{a) must converge to s as a tends to zero, the smallest solution in [0, tt] of 1 + 3cos(2s) = 0. ■ 



7 Appendix 

The following version of the inverse function theorem is used in the argument of Proposition ^. 

Theorem 3 Let $ : M" ^ M" be a map and K. C M" a compact set such that ^^jq : K $(/C) is bijective 
and the differential D^{x) is invertible for x G K. Then, there exists an open neighborhood U D K, such that 
is a diffeomorphism. 

Let now {'^k)k>i be a sequence of maps converging in the Cl^^ sense to $. Then, for every open set U 
with closure included in U , there exists k such that, for every k > k, ^k^^ is a diffeomorphism and, for 

every compact subset K. ofU, $(^) C $fc(f/) and lim <^^'^{v) = uniformly with respect to v G ^(^)- 

k — >oo 

Proof. Let us define, for k > 0, the following open neighborhoods of K, 

A:={xe R" : det L>$(x) ^ 0} Ak := U^eicB(x, ^ A . 

In view of the inverse function theorem, in order to conclude the proof of the first part, it is enough to show 
that for k large enough the restriction ^\Ak is one-to-one. 
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Figure 10: Comparison between the optimal synthesis for the hnear pendulum and the optimal synthesis on 
the bottom of the sphere in the case r{a) = 



We argue by contradiction. Let Xk ^ yk & Ak such that <I>(xfc) = ^{yk) Vfc. Then, up to extractions of 
subsequences, we can assume that the two sequences converge to x and y respectively. Since x,y € HkAk = JC 
and $(x) = we deduce that x = y. However, since detD^{x) ^ 0, we have that <& is bijective in a 

neighborhood of x, which contradicts the assumption ^(xfe) = ^(yfc) for fc large enough. 

The proof of the second part is similar. First, fix a subset U of \J . By the uniform convergence of D^fe to 
1^0 on every compact subset of [/, we get det Z?$fe(a;) 7^ for every x €U and k large enough. We also obtain 
that <i>fe is one-to-one with the same argument as above. For the remaining results to establish, they simply 
follow from the uniform convergence of <i?fe to $ on every compact subset of U . I 
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